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Abstract 

The origin of the knee in the energy spectrum of cosmic rays is an outstanding problem in 
astroparticle physics. Numerous mechanisms have been proposed to explain the structure 
in the all-particle spectrum. In the article basic ideas of several models are summarized, 
including diffusive acceleration of cosmic rays in shock fronts, acceleration via cannon- 
balls, leakage from the Galaxy, interactions with background particles in the interstellar 
medium, as well as new high-energy interactions in the atmosphere. The calculated en- 
ergy spectra and mean logarithmic masses are compiled and compared to results from 
direct and indirect measurements. 
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1 Introduction 

Cosmic-ray particles — fully ionized atomic nu- 
clei — are the only matter which reaches the 
earth from outside the solar system and are ac- 
cessible for direct investigations. The energy 
of these particles can exceed even lO^o eV, but 
their acceleration mechanisms and the sites of 
origin are still under discussion. The energy 
spectrum is of non-thermal origin and follows 
a power law dN/dE cx E'^ over many orders of 
magnitude. The spectrum steepens at energies 
around 4 PeV from a spectral index 7 ~ —2.7 to 
7 ~ —3.1. This feature is commonly called the 
knee and its explanation is generally believed to 
be a corner stone in understanding the origin of 
cosmic rays, providing answers to one of the key 
questions of astroparticle physics. 



Many approaches have been discussed in the 
literature to describe the origin, acceleration, 
and propagation of cosmic rays and several 
mechanisms have been proposed to explain the 
knee. In this article the basic ideas of theoretical 
models are briefly mentioned and the resulting 
energy spectra and mass compositions are de- 
scribed in section 01 

The steeply falling cosmic-ray flux as func- 
tion of energy requires increasing apertures 
and exposure times. Presently, above 1 PeV 
cosmic rays are experimentally accessible in 
ground-based installations only. These detec- 
tors do not measure the primary particles, in- 
stead secondary particles are measured, pro- 
duced in high-energy interactions in the atmo- 
sphere, forming an extensive air shower. This 
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makes the interpretation of the indirect mea- 
surements very difficult and the results obtained 
depend on the understanding of high-energy in- 
teractions in the atmosphere. The actual exper- 
imental status is briefly discussed in section |21 

2 Experimental results and the 

poly-gonato model 

Recently, the results of direct measurements of 
energy spectra of individual elements and indi- 
rect measurements of the all-particle spectrum 
of cosmic rays have been investigated by the au- 
thor PQ. It has been shown, that the observed 
energy spectrum and mass composition can be 
described consistently in the energy range from 
10 GeV up to at least 100 PeV, adopting a phe- 
nomenological model, the poly-gonato model. 

The energy spectra of individual nuclei are as- 
sumed to follow power laws with a cut-off at high 
energies. For the cut-off different approaches 
have been tried, the cut-off energy has been 
taken proportional to either the nuclear charge 
or the nuclear mass of the individual elements. 
In addition, also a constant cut-off energy has 
been assumed. Two possibilities for the spec- 
tral shape above the cut-off energy have been 
investigated, a common spectral slope 7c for all 
elements and a constant difference A7 of the 
spectral index before and above the cut-off en- 
ergy. 

Individual spectra of elements, as directly 
measured at the top of the atmosphere for ener- 
gies below 1 PeV, are extrapolated to high en- 
ergies and fitted to results from air-shower mea- 
surements. For these indirect measurements the 
individual energy spectra had been normalized 
to match the all-particle spectrum of direct mea- 
surements at 1 PeV. It turned out, that only 
small energy shifts were necessary, all within the 
quoted errors — but result all air-shower 
measurements yield consistent all-particle en- 
ergy spectra. 

It turned out that a cut-off proportional to 
the nuclear charge describes the data best. The 
modeling of the cut-off behavior seems not to be 
crucial for the description of the measurements, 
the two assumptions of a common 7c and a com- 
mon A7 yield almost the same results. The en- 



ergy spectra obtained are hardly distinguishable 
from each other. In addition, a significant con- 
tribution of ultra-heavy elements in the energy 
region around 10® GeV has been found. With 
the poly-gonato model two prominent features 
in the cosmic-ray spectrum can be explained. 
The knee at 4.5 PeV follows from subsequent 
cut-offs for individual elements, starting with 
protons, and the second knee at 400 PeV is due 
to the end of the galactic component. 

In the past air shower experiments mostly de- 
rived all-particle energy spectra from their mea- 
surements. However, recently also results on the 
energy spectra for groups of elements in the PeV 
region have been derived from indirect measure- 
ments. 

The KASCADE group performed systematic 
studies to evaluate the infiuence of different 
hadronic interaction models used in the simu- 
lations to interpret the data on the resulting 
spectra for elemental groups |21 01 IH El El • Two 
sets of preliminary spectra, derived from the ob- 
servation of the electromagnetic and muonic air 
shower components, using CORSIKA _7j with 
the hadronic interaction models QGSJET and 
SIBYLL are compiled in Figure ^ In order to 
compare the results with other measurements 
and theoretical predictions, the flux of the all- 
particle spectrum has been normalized to the 
average measured flux, as obtained by many ex- 
periments , represented by the solid line in the 
flgure. For the normalization, the energy scale 
has been slightly adjusted by 6e = —15% for 
both interaction models used. For comparison, 
also spectra according to the poly-gonato model 
are given. 

As can be seen in the figure, the flux for ele- 
mental groups depends on the model used. The 
KASCADE group emphasizes that at present 
there are systematic differences between mea- 
sured and simulated observables, which result 
in the ambiguities of the spectra for elemental 
groups. These conclusions apply similarly also 
to other experiments. A correct deconvolution 
of energy spectra for elemental groups requires 
a precise knowledge of the hadronic interactions 
in the atmosphere. It seems that the interac- 
tion models presently used do not describe the 
measurements with the high precision required. 
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Figure 1: Preliminary energy spectra for groups of elements as obtained by the KASCADE experiment 
investigating the electromagnetic and muonic components [Hj, as well as single hadrons [Sj, using 
CORSIKA simulations with the hadronic interaction models QGSJET (left-hand side) and SIBYLL 
(right-hand side) to interpret the data. The energy scale of the data has been readjusted, see text. 
The all-particle flux is represented by the filled symbols. The solid line gives the average all-particle 
flux from several experiments The dashed and dotted lines indicate spectra according to the 
poly-gonato model for elemental groups with nuclear charge numbers as listed. 



Nevertheless, analyzing the flux for individ- 
ual elemental groups, both sets of spectra in- 
dicate a rigidity dependent cut-off. While the 
absolute flux values obtained are different com- 
pared to the poly-gonato model, the cut-off be- 
havior for the individual groups is very similar 
to the spectra according to the phenomenologi- 
cal approach. Also presented is the flux of pro- 
tons obtained from an analysis of unaccompa- 
nied hadrons |H], which is found to be compat- 
ible with the extrapolations of the direct mea- 
surements below the knee. 

The results of the EAS-TOP experiment are 
summarized in Figure |2 The proton spec- 
trum as derived from unaccompanied hadrons 
[S] agrees well with the extrapolation of the di- 
rect measurements. The combined analysis of 
EAS-TOP Cerenkov and MACRO muon data 
[Tfl] have been normalized {Se = —15%) to 
the proton+helium flux of direct measurements, 
otherwise the flux for p+He-|-CNO would exceed 
the average measured all-particle spectrum. The 
results obtained are compatible with the extrap- 
olation of the direct measurements. The spectra 
derived from the number of electrons and muons 
(scaled hy 5e = —5%) are given for two cases 



im I12j. In the first, the light component con- 
sists of protons only and in the second scenario it 
is a mixture of protons and helium nuclei. The 
range obtained is indicated by the two sets of 
points for each mass group. 

The HEGRA collaboration derived a sum- 
spectrum for the proton and helium compo- 
nent as well as an all-particle spectrum JH]) as 
shown in FigureOl The all-particle spectrum has 
been normalized to the average measured flux 
{5e = —10%). Below 1 PeV the flux for protons 
and helium is lower than the extrapolations of 
the direct measurements, while at higher ener- 
gies the data compare well to the poly-gonato 
model. 

The figure also reviews the spectra obtained 
by the Tibet group HH UHl d , the all-particle 
spectrum has been normalized to the average 
measured spectrum {5e = —10%). The spec- 
tra for protons and helium nuclei are below the 
extrapolations of the direct measurements using 
single power laws. However, one has to keep in 
mind that these results are obtained with a very 
limited number of events observed. As exam- 
ple, the results of 2003 ^Hl are based on less 
than 200 registered events. On the other hand. 
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Figure 2: Energy spectra for groups of elements 
as obtained by the EAS-TOP experiment from 
the analysis of single hadrons the electro- 
magnetic and muonic components ^DE]) and 
a combined analysis of Cerenkov light and high- 
energy muons registered by the MACRO detec- 
tor . The energy scale of the data has been 
readjusted, see text. The all-particle flux is rep- 
resented by the filled symbols, see also caption 
of Figure ^ 



the experiments mentioned previously use data 
with high statistics, e.g. the KASCADE group 
uses O(IO^) events for their analysis. 

It is worth mentioning that all energy shifts 
applied have a negative sign, i.e. the energy is 
overestimated by the hadronic interaction mod- 
els used to interpret the measured data. A pos- 
sible explanation for this behavior could be a 
too large increase of the inelastic cross sections 
as function of energy, as pointed out in [T7j . 

It should be emphasized that the poly-gonato 
model is based on individual element spectra 
from direct measurements and the all-particle 
spectrum from indirect measurements only. The 
energy spectra for elemental groups derived from 
air shower measurements as just discussed are 
not taken into account. Hence, the poly-gonato 
model can be treated as independent result and 
it is interesting to realize that the predicted 
spectra are compatible with the deconvoluted 
spectra from indirect measurements. 

However, the present experimental status con- 
cerning energy spectra for elemental groups in 
the PeV region is not yet conclusive. Additional 
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Figure 3: Energy spectra for groups of ele- 
ments as obtained by the HEGRA jJSj and Tibet 
^JEIE] groups. The energy scale of the data 
has been readjusted, see text. The solid line 
gives the average all-particle flux from several 
experiments p. The dashed and dotted lines 
indicate spectra according to the poly-gonato 
model for protons and helium nuclei. 



information from a robust observable like the 
mean logarithmic mass would be very helpful 
to evaluate the predictions of various theoret- 
ical models of the knee. The mean logarith- 
mic mass, defined as {In A) = ^rjln^j, with 
the relative fraction of the element with mass 
number Ai is an often used quantity to charac- 
terize the cosmic-ray mass composition. 

In the earlier article |1| experimental (In^) 
values have been investigated. A systematic dif- 
ference has been found in the mean logarithmic 
mass calculated from experiments measuring 
particle distributions at ground level and detec- 
tors observing the average depth of the shower 
maximum. A recent investigation showed, that 
the discrepancy can be reduced, if the logarith- 
mic increase of the inelastic cross-sections in the 
model QGSJET is lowered and the elasticity of 
the interactions is increased [T7| . 

Taking the experimental values from ex- 
periments measuring particle distributions at 
ground level ^ and results from experiments 
observing the average depth of the shower maxi- 
mum according to an interpretation with a mod- 
ification of QGSJET (model 3a in ^Zj), the aver- 
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Figure 4: Mean logarithmic mass derived from many experiments (shaded area) compared to {In A) 
as obtained with different models. Upper left panel: Acceleration in supernova remnants as described 
by Berezhko and Ksenofontov jlH], Kobayakawa et al. jlOj, Stanev et al. and Sveshnikova [2T] . 
Upper right panel: Source and acceleration related models by Eriykin and Wolfendale [Hj, Plaga 
as well as Volk and Zirakashvili f23j. Lower left panel: Diffusion models by Kalmykov et al. Il2l, Ogio 
and Kakimoto ^28j, as well as Roulet et al. |29j . Lower right panel: Propagation models by Dova et 
al. [HJ, Lagutin et al. [2^], Swordy [23, and Tkaczyk . In addition, the range of (In^) for the 
extrapolation of direct measurements according to the poly-gonato model is indicated as dotted grey 
lines, see text. 



age mean logarithmic mass has been calculated. 
For details on the experiments and references 
the reader is referred to El • The mean value 
and the variation as function of primary energy 
are given in Tabled for reference. The average 
experimental values are also presented graphi- 
cally as shaded area in Figure |3 together with 
the predictions of various models as described 
below in section |31 

The values calculated with the poly-gonato 



model are presented in Figure 21 in addition as 
grey lines. The inclusion of ultra- heavy elements 
has only a minor effect on (In^), at 10^ GeV the 
values with and without ultra-heavy elements 
differ only by A(lnA) « 0.2. The dotted lines 
depict the systematic errors of the model re- 
sulting from the uncertainties of the direct mea- 
surements and the errors of the fit parameters 
of the poly-gonato model. In order to include 
also systematic effects caused by the assump- 
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Table 1: Average mean logarithmic mass (In^) 
as obtained from many experiments as function 
of primary energy £^0 [GeV] . fi^in denotes the 
variation, not the error of the mean value. 
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tions made in the model, {in A) has been cal- 
culated for 1 < Z < 28 and 1 < Z < 92, 
with a common 7c and a common A7, as well 
as with and without ad-hoc component. The 
largest and smallest values obtained with these 
combinations are shown in the figure, represent- 
ing the total systematic error. The errors are 
asymmetric, since compared to heavy elements 
the energy spectra of light elements are specified 
more precisely by direct measurement. 

It should be mentioned that at energies be- 
low 1 PeV the indirect observations obtain a 
lighter mass composition than the direct mea- 
surements. The (hi A) range of the latter is 
bounded by the dotted grey lines. 

3 Theoretical models for the 

knee 

Proposed explanations in the literature for the 
knee may be divided into four categories. The 
first three discuss astrophysical reasons, argu- 
ing that the knee is an intrinsic property of the 
energy spectrum, while in the last class the au- 
thors consider new particle physics processes in 
the atmosphere. In these theories the knee does 



not exist in the primary energy spectrum, but is 
only an effect of observing extensive air showers 
in the atmosphere: 

1) Models in the first category relate the 
knee to the acceleration process. The stan- 
dard approach of particle acceleration in shock 
fronts from supernova explosions and sev- 
eral extensions are discussed in the models 
by Berezhko and Ksenofontov Stanev et 
al. ^Hl, Kobayakawa et al. [10], and Sveshnikova 
|21j . The maximum energy attained is taken to 
be responsible for the knee. Also source related 
is the proposal of Erlykin and Wolfendale ^22] . 
namely that a nearby supernova as a single 
source of cosmic rays causes the structure in 
the energy spectrum. Volk and Zirakashvili 
consider reacceleration by spiral shocks in the 
galactic wind. Plaga |21] proposes cosmic rays 
being accelerated by cannon balls ejected into 
the galactic halo. 

2) Models of the second category connect 
the knee with leakage of cosmic rays from the 
Galaxy. Swordy combines diffusive shock 
acceleration with an energy dependent propaga- 
tion pathlength in a Leaky Box model. In the 
models by Lagutin et al. j^H] Ptuskin et al. j27j , 
Ogio and Kakimoto |2H]) as well as Roulet et al. 
j29j the knee follows from diffusive propagation 
of cosmic-ray particles through the Galaxy. 

3) Interactions of cosmic rays with back- 
ground particles in the Galaxy are considered 
as origin for the Jcnee in the third category. 
Tkaczyk [Hni suggests diffusive propagation in 
combination with photo-disintegration. Interac- 
tions of cosmic rays with background photons 
are proposed by Dova et al. Candia et 
al. discuss the interaction of cosmic rays 
with the neutrino background. 

4) A last class of theories accounts the air 
shower development in the atmosphere for the 
knee. The basic idea is that a new type of in- 
teraction transfers energy to a component not 
observable (or not yet observed) in air shower 
experiments. The threshold for these new inter- 
actions is in the knee region, above the energy 
of todays collider experiments. In the model 
by Kazanas and Nicolaidis the energy is trans- 
fered into techni-hadrons, lightest supersymmet- 
ric particles [HH]) and gravitons |3l] . 
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In the original publications different represen- 
tations for the theoretical spectra are used, in- 
cluding different sets of mass groups and en- 
ergy multipliers for the flux values. In order 
to provide an overview on the predicted spec- 
tra and to allow an easy comparison, all re- 
sults are presented in a standard format with 
the ordinate being multiplied by E'^'^. The en- 
ergy spectra are given for the following five mass 
groups: protons (Z=l), helium (Z=2), interme- 
diate: lithium - fluorine (Z=3~9), heavy: neon - 
titanium (Z=10-22), and very heavy: vanadium 
- uranium (Z=23-92). 

To provide some orientation for the compar- 
ison of different models, in the figures depict- 
ing the results of the individual models also cal- 
culations according to the poly-gonato model 
are shown. As has been discussed in section |21 
these spectra are reasonable approximations of 
the spectra for elemental groups as derived from 
air shower measurements. 

3.1 Acceleration in supernova rem- 
nants 

The first four models deal with different varia- 
tions of diffusive shock acceleration of cosmic 
rays in shock fronts. The maximum energy 
reached in this process is related to the knee 
in the energy spectrum. The calculations by 
Berezhko and Ksenofontov TH* are based on the 
nonlinear kinetic theory of cosmic-ray accelera- 
tion in supernova remnants. The mechanical en- 
ergy released in a supernova explosion is found 
in the kinetic energy of the expanding shell of 
ejected matter. The cosmic-ray acceleration is a 
very efficient process and more than 20% of this 
energy is transfered to ionized particles. The ac- 
celeration is achieved in the following way: The 
diffusive propagation of high-energy particles al- 
lows them to cross the shock front many times 
and every pair of consecutive crossings of the 
shock front increases the energy of the particles. 
In linear approximation this process generates 
a power-law momentum spectrum at a planar 
shock front. Due to the high acceleration ef- 
ficiency and the hardness of the spectrum the 
structure of the shock wave is modified by the re- 
ciprocal influence of cosmic rays on the medium 
and, in consequence, the cosmic-ray spectrum is 



altered. 

The particles need some minimum velocity in 
order to cross the shock front. This velocity de- 
termines the injection rate of particles into the 
acceleration regime. Superthermal particles can 
efficiently be injected and accelerated with high 
efficiency. The injection efficiency is expected to 
be an increasing function of the mass to charge 
ratio (A/Z) of the nucleus considered, i.e. heavy 
elements are accelerated more efficiently and ob- 
tain a harder spectrum. The energy limit of cos- 
mic rays accelerated in supernova remnants is 
controlled by geometrical factors of the expand- 
ing shock front. Taking into account preacceler- 
ation in the wind of the predecessor star maxi- 
mum energies up to Z • 10^^ eV can be achieved. 
In regions modified by the presupernova wind 
the magnetic fields are larger than in the unper- 
turbed interstellar medium. For the calculations 
the density of elements heavier than hydrogen 
in the interstellar medium is assumed to be pro- 
portional to the relative abundance of elements 
in the local, solar region of the Galaxy. The 
expansion of the shock wave into a hot inter- 
stellar medium, which is assumed to be signifi- 
cantly modified by presupernova winds of type 
lb and type II supernovae, with a magnetic field 
Bq = 12 has been taken into account with 
an injection rate t] = 5 ■ 10~^. 

The resulting energy spectra for five elemen- 
tal groups are plotted in Figure \^ (symbols). 
The knee in the all-particle spectrum is due to 
the charge dependence of the maximum energy 
achieved in the acceleration process. The proton 
component exhibits a knee at 1 PeV with a rel- 
atively small change of the spectral index which 
keeps the protons to be the dominant compo- 
nent up to 10^ GeV. Elemental spectra accord- 
ing to the poly-gonato model are represented 
by the lines. They show a steeper decrease of 
the flux above the knee for all elements. Below 
1 PeV the spectra for He and the CNO-group 
of the model by Berezhko &: Ksenofontov agree 
with the poly-gonato model, but the heavy el- 
ements are suppressed in favor of the protons. 
On the other hand, the all-particle spectra (solid 
line and filled dots) are well compatible in shape 
and absolute normalization. 

The mean logarithmic mass has been calcu- 
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Figure 5: Energy spectra for five groups of elements according to tlie model by Berezhko & Ksenofontov 
|18j (open symbols) and the poly-gonato model (dashed and dotted lines), the numbers indicate the 
range of nuclear charge for each group. Filled points and the full line represent the all-particle spectra. 



lated using the five mass groups, the result is 
shown in Figure \^ as function of energy^ . The 
values exhibit a slight decrease up to about 
10^ GeV to a minimum value around (In A) = 
1.2. At higher energies {In A) increases and 
reaches (In^) = 2.1 at 10^ GeV. The rise of 
(In^) beyond the knee is relatively modest for 
the model in question, mainly caused by the 
moderate flux decrease beyond the individual 
knees. In the whole energy range the (In A) val- 
ues are at the lower edge or even below the range 
of measured values. 

3.2 Acceleration by supernova shocks 

A threefold origin of energetic cosmic rays is pro- 
posed by Stanev, Biermann, and Gaisser jT^. In 
their model particles are accelerated at three dif- 
ferent main sites j35j : 

1) The explosions of normal supernovae into 
an approximately homogeneous interstellar 
medium drive blast waves which can accelerate 
protons to about 10^ GeV. Particles are accel- 
erated continously during the expansion of the 
spherical shock-wave, with the highest particle 

^Only (A) versus energy is given in the original article. 



energy reached at the end of the Sedov phase. 

2) Explosions of stars into their former stellar 
wind, like that of Wolf Rayet stars, accelerate 
particles to higher energies. The maximum en- 
ergy attained depends linearly on the magnetic 
field and maximum energies Emax = 9 • lO"^ GeV 
for protons and E^ax = 3 • 10^ GeV for iron nu- 
clei are reached. 

3) For energies exceeding 10*^ GeV an extra- 
galactic component is introduced by the au- 
thors. The hot spots of Fanaroff Riley class II 
radio galaxies are assumed to produce particles 
with energies up to 10^^ GeV. 

In this model the main fraction of galactic 
cosmic rays above about 10 TeV are assumed 
to be accelerated by shocks that travel down a 
steady stellar wind with a Parker spiral struc- 
ture for the magnetic field. Like in other first 
order Fermi theories the basic idea for the en- 
ergy gain is the cyclic crossing of the shock wave 
from upstream to downstream and back. In 
most regions the shock normal is perpendicu- 
lar to the magnetic field. But in small regions 
around the poles (about 1% of Att sr) the direc- 
tion of the propagation of the shock is parallel 
to the magnetic field which yields harder spectra 




Figure 6: Energy spectra for five groups of elements according to the model by Stanev, Biermann, 
and Gaisser JHI (open symbols) and the poly-gonato model (dashed and dotted lines). The numbers 
indicate the range of nuclear charge for the elemental groups. Filled points show the sum spectrum 
for all groups according to the theoretical calculations. The sum spectrum for all elements according 
to the poly-gonato model is represented by a solid line, above 10^ GeV a dashed line represents the 
average experimental flux spectrum (see j]P). 



for the accelerated particles. The combination of 
the polar cap with the rest of the stellar hemi- 
sphere might lead to a situation where up to 
about 10^ GeV the entire hemisphere excluding 
the polar cap dominates, while from 10^ GeV up 
to the icnee the polar cap begins to contribute 
appreciably. At knee energies the polar cap com- 
ponent contributes equally to the remainder of 
the 47r sr solid angle. This fact is essentially to 
describe a sharp bend in the all-particle spec- 
trum as it has been observed by some experi- 
ments, e.g. Akeno ISHl- The logical consequence 
of the model is that the chemical composition at 
the knee changes in a way, that the gyroradii of 
the particles at the spectral break are identical. 
This implies that the fluxes of different nuclei 
break off according to their charge Z. 

Taking the chemical composition as measured 
directly by experiments up to 100 TeV ener- 
gies, Stanev et al. fit the shower size spectra 
obtained by the Akeno experiment for different 
zenith angles and derive primary energy spec- 
tra for groups of elements. The spectra ob- 



tained for five mass groups are reproduced in 
Figure El The proton flux starts to deviate from 
a power law at 10^ GeV caused by the energy 
limit reached for the Sedov phase expansions in 
the interstellar medium and the flux decreases 
up to 10^ GeV. The increase of the proton flux 
above this energy is caused by the extragalactic 
component introduced. Elements heavier than 
hydrogen exhibit a sharp, rigidity dependent 
bend. 

The relative distinct structure of the knee is 
due to the contribution of the polar cap com- 
ponent. Beyond the knee the calculated spectra 
are flatter as compared to the poly-gonato model 
and the contribution of light elements is larger. 
The sharp knees for the individual element spec- 
tra result in a pronounced structure in the all- 
particle spectrum. The model has been tuned 
to describe the sharp knee seen by early exper- 
iments like MSU and Akeno. Recent measure- 
ments yield a smoother transition at the knee, 
see e.g. Figure 15 in 1. A fit to the more recent 
data would probably result in a smaller contri- 
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Figure 7: Energy spectra for groups of elements for the model by Kobayakawa et al. j2L)j (symbols) 
and the poly-gonato model (lines), see caption Figure El 



bution of the polar cap component. 

The mean logarithmic mass calculated from 
the individual element spectra, shown in Fig- 
ure lU exhibits an increase up to (lnj4) ~ 3.1 at 
an energy Eq ~ 10'' GeV. The early increase is 
caused by the early decrease of the proton flux 
mentioned. The (In^) values predicted by the 
model increase faster as function of energy than 
the observations, but reach their maximum at 
lower energies and smaller (In A) values are pre- 
dicted at the maximum. The first peak is due to 
the Jciiees in the spectra for the heavy elements 
(Ne and Fe groups) at energies around 10'' GeV. 
The second peak and the decreasing (lnj4) val- 
ues are caused by the extragalactic proton com- 
ponent. 

3.3 Acceleration by oblique shocks 

A slightly modified version of the diffusive ac- 
celeration of particles in supernova remnants is 
considered by Kobayakawa, Sato, and Samura 
j2()j . Standard first order Fermi acceleration in 
supernova remnants — with the shock normal 
being perpendicular to the magnetic field lines 
— is extended for magnetic fields with arbitrary 
angles to the velocity of the shock front. The ba- 



sic idea is that particles are accelerated to larger 
energies in oblique shocks as compared to par- 
allel shocks. 

The ejected material from a supernova explo- 
sion expands into the interstellar medium driv- 
ing a shock wave. The shock acceleration may 
be effective until the expanding shell sweeps up 
its own mass of the interstellar medium. The 
ejected mass (~ IOMq) with large momenta 
moves out freely at constant velocity during the 
free expansion phase which is assumed to be in 
the order of a few hundred years for ejecta ex- 
panding at a velocity ~ 10^ cm/s into a medium 
of an average density ~ 1 proton/cm'^. The 
shape of the shock front is supposed to be al- 
most spherically symmetric. The directions of 
the interstellar magnetic field lines will be, over 
a wide range, nearly at random rather than well 
aligned. Therefore, it is assumed that the field 
lines meet the shock front at random angles and 
the cosines of the angles are distributed uni- 
formly. 

The injection efficiency into the acceleration 
regime is investigated as function of the angle 
between the magnetic fields and the normal of 
the shock front. It is found that the injection ef- 
ficiency decreases drastically as the shock obliq- 
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uity increases. The simple relation for the ef- 
ficiency e{ri) = 7] is assumed, rj being the co- 
sine of the angle between the magnetic field and 
the expansion direction of the shock front. On 
the other hand, the obliquity makes the energy 
spectrum harder, particles are accelerated more 
efficiently. 

For oblique shocks, reflection at the shock 
front is more important than diffusion in down- 
stream regions, so that a more rapid acceler- 
ation can occur than at parallel shocks. The 
maximum energies of accelerated particles can 
be two or even three orders of magnitude larger 
in the case of a quasi-perpendicular shock than 
in a parallel one. The maximum energy per par- 
ticle reached during the acceleration process is 
mainly limited by the finite lifetime of the su- 
pernova blast wave. In the article a rigidity de- 
pendent maximum energy 

Emax oc Z • 2.5 • 10^6 eV 

_B . u . fi . (1) 

is derived, with the speed of the shock wave n, 
the upstream magnetic field B — being about 
ten times larger than conventionally estimated 
— , and /(??) being an angle dependent acceler- 
ation efficiency. For parallel, strong shocks {rj = 
1) a maximum energy Ecrit = Z ■ 1.25 • 10"*^^ eV is 
obtained. This value is exceeded by more than a 
factor of 1400 for shocks with extreme obliquity. 

The authors fit experimental data for several 
groups of elements. The spectral indices of the 
power laws as well as the absolute flux are ad- 
justed at 1 TeV. The resulting spectra are ex- 
trapolated to high energies as presented in Fig- 
ure d The individual knees result from the de- 
pendence of the maximum energy reached on the 
angle between the shock and the magnetic field. 
The spectra exhibit a rigidity dependent icnee. 
A comparison with the poly-gonato model ex- 
hibits some differences, but the overall agree- 
ment is quite good. However, above the knees 
the model yields slopes less steep and overesti- 
mates the all-particle flux. 

The cut-off behavior of the iron-group is the 
reason for the strong increase of the mean loga- 
rithmic mass as function of energy, as shown in 
Figure 13 The model predictions exhibit an in- 
crease up to (Inyl) « 3.5 at 10^ GeV, very close 



to the maximum of the poly-gonato model. The 
(In A) values are at the upper border or even 
above the experimental values in the whole en- 
ergy range. 

3.4 Acceleration by a variety of su- 
pernovae 

In a recent article Sveshnikova revised the stan- 
dard approach of acceleration of cosmic rays in 
shock fronts of supernovae j21j . The new ap- 
proach to understand acceleration in supernova 
remnants (SNRs) is based on recent astronom- 
ical observations of supernovae [SZIISHI- They 
indicate that core collapse supernovae prove to 
comprise the most common general class of ex- 
ploding stars in the universe and they occur in 
a great variety of flavors. In particular, super- 
novae have been observed with unusual large ex- 
pansion velocities, indicating that these objects 
are hyper-energetic, the so called hypernovae. 

The upper limit of acceleration in SNR is de- 
termined essentially by the product of the shock 
radius Rsk, shock velocity Vsk = Vgk/WOO km/s, 
ejected mass Mgj, remnant age Tsnr, and explo- 
sion energy £^51 = Egnr/W^^ erg. All these val- 
ues are connected to each other and vary from 
explosion to explosion. The cut-off energy per 
particle Smax[TeV] can be expressed by a simple 
formula if only the Sedov Phase of SNR expan- 
sion is considered [30] 




= Z-E^,,{B,nH){E,i-Vsk)^^ 

The maximum energy depends on three factors, 
the charge Z of the nucleus, the properties of the 
interstellar medium where the SNR is expand- 
ing (strength of magnetic field B and density of 
protons uh), and on the energy of the explosion 
as well as the velocity of the shock. 

A usual way to increase E^ax is to enlarge 
the magnetic field B, see e.g. the approach 
by Kobayakawa et al. (section II13|) . On the 
other hand, as consequence of the recent obser- 
vations also the parameters of the supernova ex- 
plosion can be varied. This is the basic idea of 
Sveshnikova, to draw up a scenario in which the 
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Figure 8: Energy spectra for groups of elements for the model by Sveshnikova pij (open symbols) 
and the poly-gonato model, see caption of Figure El 



maximum energy reached in SNR acceleration 
is in the icnee-region of the cosmic-ray spectrum 
(~ 4 PeV), using only the standard model of 
cosmic-ray acceleration and the latest data on 
super novae explosions. 

Based on recent observations the distribution 
of explosion energies £^51 and their rates of oc- 
currence in the galaxy are estimated. The spec- 
trum of cosmic rays in each explosion is approx- 
imated by a composite power law 



^(E) = ^oE-'f with 
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E > Ejriax/^ 



(3) 



The observed spectrum in the Galaxy is ob- 
tained as sum over all different types of super- 
novae explosions, integrated over the distribu- 
tion of explosion energies within each supernova 
group. The mass composition of cosmic rays 
has been assumed according to direct measure- 
ments. As result energy spectra at the source 
are obtained for groups of elements. In order to 
take into account propagation effects the spectra 
given in 21] have been multiplied with Eq^'^^, 
as suggested in the article. The flux has been 
normalized to the average measured all-particle 



spectrum at low energies and the resulting spec- 
tra are depicted in Figure El 

Since the spectra for elemental groups are the 
result of many supernovae explosions, they can 
not be described by a simple function with two 
power laws. Instead, they exhibit a more com- 
plicated structure, the cut-off above the indi- 
vidual knees is characterized by two subsequent 
steps, caused by contributions of individual su- 
pernovae classes. The the first step for each 
elemental group occurs at about the same en- 
ergy as in the poly-gonato model. At energies 
larger than the individual cut-offs, the spectra 
predicted by Sveshnikova are significantly flatter 
due to the contribution of very energetic SNRs. 

The knee in the all-particle spectrum is caused 
by the cut-off of the proton component. The 
all-particle spectrum above the knee is formed 
by the sum of coupled steps. The change of 
the slope beyond the ioiee (in the interval from 
3 PeV to Zpe ■ 3 PeV) is determined by the frac- 
tion of iron nuclei in the chemical composition 
of cosmic rays before the knee. The maximal 
energy of Galactic cosmic rays is determined by 
iron nuclei generated in the most energetic su- 
pernovae (type SNIIn) and corresponds to about 
10^ GeV. 
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Figure 9: Energy spectra for groups of elements for the model by Eriykin & Wolfendale j22j (open 
symbols) and the poly-gonato model, see caption of Figure El 



The iron component is clearly dominant at 
high energies, which leads to a strong increase 
of the mean logarithmic mass, as can be in- 
ferred from Figure |1J A slight structure is visible 
around 10 PeV, at this energy the light nuclei 
(protons and helium) are replaced by heavy nu- 
clei as the major component in cosmic rays. The 
mean logarithmic mass increases to values above 
{In A) = 3.5 at 10^ GeV. Additionally, a vari- 
ant is discussed, where cosmic rays from SNIIn 
are enriched by heavy nuclei. This results in 
a heavier mass composition. In the region be- 
tween 10^ and 10® GeV (In A) is increased up 
to 1 unit. The model predictions for both vari- 
ants are well compatible with the average ex- 
perimental values. The latter scenario exhibits 
good agreement with the poly-gonato model. 

3.5 The single-source model 

A different approach is used by Eriykin and 
Wolfendale 22 . In addition to a background 
caused by many undefined sources the authors 
presume a single supernova remnant from a re- 
cent and nearby explosion as additional source 
of cosmic rays. The knee in the all-particle 
spectrum is supposed to have a two-kink struc- 



ture related to the cut-offs of oxygen and iron 
nuclei from the single source. The authors 
use shower size spectra of the electromagnetic, 
muonic, hadronic and Gerenkov components of 
extensive air showers. The shower size spec- 
tra are normalized to the same knee position. 
The proposed twofold structure is not visible in 
all individual shower size spectra considered, it 
shows up after the normalization and rebinning 
procedure, combining results from several exper- 
iments. 

In the publication energy spectra for groups 
of elements are given for both, the background 
and the single SNR separately. The sum of both 
components is shown in Figure 1^1 for individual 
element groups. Salient kinks in the spectra 
of the GNO and Fe groups caused by the sin- 
gle SNR produce a twofold structure in the all- 
particle spectrum at 3 • 10^ GeV and 10^ GeV. 
Such a structure is not observable in the aver- 
aged all-particle spectrum of the measurements. 
In the single-source model the knees of the el- 
ements heavier than helium are much sharper 
than in the poly-gonato model. 

The mean logarithmic mass taken from Er- 
iykin and Wolfendale is shown in Figure^ The 
early decrease of the proton component results 
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in a strong increase of (In A) at the upper border 
of the observed values. Within the present ex- 
perimental resolution, the structure at 10^ GeV 
can not be resolved. 

3.6 Recceleration in the galactic wind 

Reacceleration of cosmic-ray particles in the 
Galactic wind is discussed by Volk and Zi- 
rakashvili The wind is mainly driven by 

cosmic rays and hot gas generated in the disk. It 
reaches supersonic speeds at about 20 kpc above 
the disk, and is assumed to be very extended 
(several 100 kpc) before it ends in a termination 
shock. Due to Galactic rotation the differences 
in flow speed will lead to strong internal wind 
compressions, bounded by smooth cosmic-ray 
shocks. These shocks are assumed to reaccel- 
erate the most energetic particles from the disk 
by about two orders of magnitude in rigidity, 
ensuring a continuation of the energy spectrum 
beyond the knee up to the ankle. A fraction 
of the reaccelerated particles will return to the 
disk, filling a region around the Galactic plane 
(several tens of kpc thick) rather uniformly and 
isotropically. 

A maximum energy Emax ~ Z • 10^"^ eV is ob- 
tained and the authors conclude that the knee in 
the all-particle spectrum cannot be the result of 
the propagation process, instead it is supposed 
to be a feature of the source spectrum itself. It is 
pointed out that it is possible to explain the con- 
tinuation of the cosmic-ray spectrum above the 
knee up to the ankle in a natural way, by consid- 
ering the dynamics of the interstellar medium of 
the Galaxy and its selfconsistent extension into 
a large-scale halo by the Galactic wind. The au- 
thors conclude further that within this picture 
there is no way to produce higher energy cosmic 
rays, their sources must be of a different nature. 

For the calculations the chemical composition 
of cosmic rays has been adjusted at 9-10^'' eV to 
KASCADE measurements. The energy spectra 
derived for four elemental groups as well as the 
all-particle spectrum are depicted in Figure 
For the light elements (protons and helium) a 
icnee-like structure is clearly visible. For these 
elemental groups the decrease at high energies 
is significantly attenuated and the spectra con- 
tinue at high energies with a spectral slope sim- 



ilar to the values below the knee. This behav- 
ior differs from the poly-gonato model with its 
assumptions on the cut-off behavior, while the 
all-particle spectra agree well with each other. 

The mean logarithmic mass deduced is de- 
picted in Figure m In the model iron nuclei are 
only about a factor of two more abundant than 
protons at high energies (£^0 ~ 10^ eV). Con- 
sequently, only a modest increase of the mean 
logarithmic mass is obtained with a maximum 
value {In A)max ~ 2.6. The predictions are well 
compatible with the average experimental (In A) 
values below 10^ GeV, but may be slightly too 
low at the highest energies. 

3.7 The cannonball model 

Within the framework of the cannonball model, 
originally proposed by Dar and De Riijula flOl to 
explain gamma-ray bursts, Plaga discusses 
a mechanism for the acceleration of cosmic- 
ray hadrons. The author investigates if masses 
of baryonic plasma ("cannonballs"), ejected in 
bipolar supernova explosions, could be the uni- 
versal sources of hadronic galactic cosmic rays. 
It is assumed, that in each of the symmetric 
cannonball jets a total energy of 10^^ erg is 
ejected. Supernova explosions which eject can- 
nonballs should occur in the Galaxy at rates of 
1/(50 a) to 1/(600 a). Less than 1/40 of the 
total cannonball-energy has to be converted to 
the energy of the cosmic-ray particles in order 
to match the observed cosmic-ray luminosity. 

Two scenarios for the acceleration are 
sketched. The first is based on ultra-relativistic 
shocks in the interstellar medium and could ac- 
celerate cosmic rays up to knee energies. The 
second is based on second order Fermi accel- 
eration inside the cannonballs. For cosmic-ray 
propagation a diffusion coefficient D oc £^~", 
with a = 0.5, is adopted. No discontinuities 
in the energy dependence of the diffusion coeffi- 
cient are assumed, thus the icnee in the energy 
spectrum is related to source properties. 

The bulk of cosmic rays up to the knee 
is taken to be accelerated at ultra-relativistic 
shocks driven by cannonballs into the interstel- 
lar medium. The accelerated particles are im- 
mediately swept up into the plasmoid and re- 
main confined there until it has slowed down to 
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Figure 10: Energy spectra for groups of elements for the model by Volk and Zirakaslivili j23j (open 
symbols) and the poly-gonato model, see caption of Figure El 
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Figure 11: Energy spectra for groups of elements for the model by Plaga 24 (open symbols) and the 
poly-gonato model, see caption of Figure El 



subrelativistic speeds in the galactic halo. The 
maximum energy 

R Z-BT 



En 



3 • 10^^ eV 



(4) 



0.02 pc fiG 50 
that can be reached in this acceleration process 



is proportional to the radius of the plasmoid R, 
to its Lorentz factor F, and to the interstellar 
magnetic field B. Z \s the nuclear charge of the 
accelerated particle. A fraction of particles con- 
tained inside the plasmoid escapes by diffusive 
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processes and determines the observed cosmic- 
ray luminosity. A power law with an spectral 
index 7 = —2.2 at the source is obtained below 
the knee. 

Even higher energies could be reached inside 
the cannonballs by further acceleration of con- 
tained hadrons in the internal turbulent plasma 
via second order Fermi mechanisms. Typical ac- 
celeration times are in the order of several 10^ 
to 10^ a and energies up to Z • 10^'^ eV could 
be reached. Again, the particles are diffusively 
released, their energy spectrum follows a power 
law with an index 7 = —2.5 at the source. 

The energy spectra for groups of elements as 
obtained by Flag presented in Figure 
The individual spectra exhibit a knee propor- 
tional to the charge. However, the change in 
spectral slope is extremely soft. The individ- 
ual groups decline in the figure almost parallel, 
without any significant change of the mass com- 
position. The obtained all-particle spectrum is 
in reasonable agreement with the average mea- 
sured flux. 

In the whole energy range shown, light parti- 
cles dominate the mass composition. Therefore, 
only a small increase of the mean logarithmic 
mass is found above the knee, as can be inferred 
from Figure^ At an energy of 8 • 10^ GeV the 
maximum value reached is still below 2. The ex- 
treme light mass composition predicted by the 
model seems not to be compatible with the ex- 
perimental values. 

3.8 The minimum-pat hlength model 

The second class of models discussed describes 
the propagation of cosmic rays through the in- 
terstellar medium. In these scenarios the knee 
is a consequence of leakage of particles from the 
Galaxy. The calculations by Swordy are 
based on a Leaky Box model for the cosmic- 
ray propagation. The spectra of particles accel- 
erated by diffusive shocks are assumed to have 
the same spectral slope 7 = —2.15 for all ele- 
ments at the source with a rigidity dependent 
cut-off at a rigidity R ~ 10^^ V. Above this 
cut-off the spectra are assumed to decrease as 
dN/dR oc R~^. The pathlength Ag for escape 
from the Galaxy declines with rigidity but has 
some minimum value. As function of R it has 



been parametrized using the expression 

Ae = Ao i^—j + K (5) 

with Ao = 6 g/cm2, 6 = 0.6, and Rq = 10 GV. 
The value for the minimum pathlength A,- = 
0.013 g/cm^ has been determined by a fit to the 
measured all-particle spectrum. 

Using this set of parameters, the author gives 
fractional abundances for five elemental groups 
as function of energy. They are normalized to 
the mean values of direct measurements in the 
TeV energy range. Taking the fractional abun- 
dances and using the average all-particle spec- 
trum as obtained from many experiments, the 
energy spectra shown in Figure are calcu- 
lated. The change of slope for the individual 
spectra is very smooth in the model by Swordy. 
Above the knee the spectra are significantly flat- 
ter as compared to the poly-gonato model. The 
cut-off with R~^ seems to be to smooth. 

The individual spectra show a slight dip just 
below the individual Jcnees. This feature causes 
the dip in the mean logarithmic mass shown in 
FigurelH The (In^) values have been calculated 
from the individual spectra shown in Figure EF- 
They reach a maximum of (In A) ~ 2.5. The 
modest increase obtained is due to the relative 
flat spectra above the individual cut-offs, which 
yields to a mass composition, which is mostly at 
the lower edge of the experimental values. 

3.9 Anomalous diffusion in the 
Galaxy 

Lagutin et al. jj^] discuss the knee structure 
as being linked to anomalous diffusion of cosmic 
rays in the magnetic fields of the Galaxy. In the 
literature several models are considered to de- 
scribe the propagation of cosmic rays through 
the interstellar medium. Frequently, a concep- 
tion is used in which the propagation process 
is determined by scattering at magnetic-field in- 
homogeneities which have small-scale character. 
They can be considered as a homogeneous Pois- 
son ensemble and the propagation is described 
by a normal diffusion model, e.g. as discussed 
by Berezinsky et al. |41j . 

^The original publication gives only values for the 
mean mass. 
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Figure 12: Energy spectra for groups of elements for the model by Swordy ,25, (open symbols) and 
the poly-gonato model, see caption of Figure El 
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Figure 13: Energy spectra for groups of elements for the model by Lagutin et al. ^.2.6,1 (open symbols) 
and the poly-gonato model, see caption of Figure 1^1 



Lagutin et al. propose the existence of mul- 
tiscale structures in the Galaxy, related to fil- 
aments, shells, and clouds which are widely 
spread in the interstellar matter. A rich vari- 
ety of structures that are created in interact- 
ing phases having different properties can be 



related to the fundamental property of turbu- 
lence. The tool of fractal geometry is used to 
characterize the interstellar medium and, corre- 
spondingly, the magnetic field strength. In such 
a fractal-like interstellar medium the cosmic-ray 
propagation can not be described by normal dif- 



17 



fusion. Instead, an anomalous diffusion model 
or superdiffusion is proposed. The anomaly re- 
sults from large free paths (Levy flights) of par- 
ticles between magnetic domains. These paths 
are distributed according to an inverse power 
law. Also rest-states of motion in a trap (Levy 
waiting times) are investigated. 

The diffusion coefficient is assumed to depend 
on the rigidity of the particles as D oc DqR^ . Dq 
describes the propagation taking into account 
both, long free paths and the probability to stay 
in magnetic traps during the propagation. The 
anomalous diffusion coefficient is evaluated us- 
ing experimental data on the anisotropy of the 
cosmic-ray flux in the energy region from 10^ 
to 10^ GeV/nucleus. In the model the spec- 
tral index as function of primary energy is de- 
rived. The differential flux of all particles is 
supposed to consist of two components. The 
first one describes the contribution of nearby 
sources {d < 1 kpc), including 16 supernova 
remnants, being responsible for the high energy 
region including the knee. The second deals with 
the bulk of observed cosmic rays with energies 
from 0.1 GeV to 10 GeV from numerous distant 
sources {d > 1 kpc). 

The authors conclude that the knee in the 
primary cosmic-ray spectrum and the observed 
distinction in the spectral index of protons and 
other nuclei can be explained by superdiffusion 
propagation. The energy dependence of parti- 
cle spectra supports the hypothesis, that nearby 
supernova bursts are the source of high-energy 
cosmic rays. 

The calculated spectra are shown in Figure !!^ 
Between 1 and 50 PeV absolute flux and shape 
of the all-particle spectrum are compatible with 
the measured spectrum, but below and above 
these energies the calculations overestimate the 
experimental all-particle flux. The energy spec- 
tra for individual elements exhibit a very smooth 
behavior in the knee region, no kink in the spec- 
tra is visible and no distinct energy for the icnee 
can be specified. At low energies the fluxes of 
the individual element spectra do not agree with 
direct measurements and their extrapolation us- 
ing power laws. 

The model is characterized by a very smooth 
change of the spectral slopes. At the highest 



energies above 10^ GeV the iron-group is only 
about a factor of two more abundant than pro- 
tons. The resulting mean logarithmic mass, pre- 
sented in Figure 13 shows a slow — in the loga- 
rithmic plot almost linear — increase as function 
of energy. At 400 PeV a value of (In^) 2.6 is 
reached. 

3.10 Drift in the global regular mag- 
netic field of the Galaxy 

The next three models described (sections 13.101 
to I3.12|) deal with variations of similar basic 
ideas. They treat the diffusion of cosmic rays 
in the Galaxy taking into account the regular 
magnetic field, a random component, and anti- 
symmetric diffusion. 

Ptuskin et al. demonstrate that a drift 
(hall diffusion) of particles in the global regular 
magnetic field of the Galaxy may be responsi- 
ble for the knee. In the galactic halo a random 
field is assumed to exist simultaneously with the 
regular magnetic field, which is mainly toroidal 
and whose value is comparable with the value 
of the field in the gas disk. Particle diffusion 
in a magnetically-active medium is described by 
a tensor. The Hall diffusion, viz. the antisym- 
metrical part of the cosmic-ray diffusion tensor, 
is not effective at small energies, but it may 
control the leakage of cosmic rays out of the 
Galaxy at energies above 3 PeV. In this way the 
icnee in the cosmic-ray power law spectrum is 
obtained without resorting to any assumptions 
concerning this feature in the source spectrum 
or changes in the energy dependence of the dif- 
fusion coefficient along the field lines. 

The global structure of the regular magnetic 
field in the Galaxy is very important to describe 
the propagation of cosmic rays since it deter- 
mines the orientation of the diffusion tensor. Re- 
liable information on the regular field in the disk 
is available, revealing that it is predominantly 
toroidal. Observations indicate that the field 
lines of the regular field are closed circles and 
the field changes sign with radial distance after 
each 3 kpc. It has the same direction above and 
below the galactic plane. 

Only scant information on the magnetic field 
beyond the galactic disk — i.e. at distances 
> 500 pc — exists. The orientation of the 
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Figure 14: Energy spectra for groups of elements for the model by Ptuskin et al. ^2 adapted to 
experimental data by Kalmykov and Pavlov ^42^ (symbols) and the poly-gonato model, see caption of 
Figure El 



regular field in the Galactic halo is not known, 
but the toroidal components seem to dominate. 
Ptuskin et al. propose two scenarios, a flat halo 
model — it is assumed that the propagation re- 
gion of cosmic rays in the Galaxy has the form 
of a cylinder and the height of the cylinder is 
much smaller than its radius — and a large halo 
model. The latter is supported by radio astro- 
nomical data which indicate that our Galaxy has 
an extended halo with a height of about 10 kpc 
comparable to the galactic radius ~ 20 kpc. Nu- 
merical calculations of the cosmic-ray propaga- 
tion in such an environment, assuming a power- 
law energy spectrum with a constant spectral 
slope at the source, yield a cosmic ray energy 
spectrum with a knee at 3 PeV. 

At energies in the GeV region the Hall dif- 
fusion is insignificant and plays a minor part 
in cosmic-ray leakage from the Galaxy. The 
Hall diffusion coefficient increases rapidly with 
energy and begins to dominate the slowly- 
increasing usual diffusion at ~ 3 PeV, thus gen- 
erating the knee. Diffusion and drift of different 
cosmic-ray nuclei at ultrarelativistic energies de- 
pend on the energy per unit charge E/Z. Thus, 
the position of the knee for each kind of nucleus 



is proportional to Z. 



Based on the large halo model Kalmykov and 
Pavlov j42j have fitted the primary mass com- 
position in order to obtain best agreement with 
the experimental size spectrum of the MSU ar- 
ray. The resulting energy spectra are shown in 
Figure ITU for five mass groups as well as the all- 
particle spectrum. The flux given by the authors 
has been normalized to the average experimen- 
tal all-particle flux at 1 PeV. A rigidity depen- 
dent cut-off for the individual element spectra is 
clearly visible. 



Except for the absolute normalization, the 
shape of all elemental groups is very similar 
for this approach and the poly-gonato model. 
Above 30 PeV the iron group is clearly domi- 
nant. This is also reflected by the pronounced 
increase of the mean logarithmic mass, as can 
be inferred from Figure [l] The model predicts a 
very strong increase of (In A) within a small en- 
ergy interval. At 10^ GeV no other model yields 
such high values of (Invl). 
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Figure 15: Energy spectra for groups of elements for the model by Ogfo and Kakimoto j28j (open 
symbols) and the poly-gonato model, see caption of Figure El 



3.11 Diffusion in turbulent galactic 
magnetic fields 

Ogio and Kakimoto |2H| consider the regular 
magnetic field in the Galaxy following the di- 
rection of the spiral arms. In addition, irreg- 
ularities of roughly the same strength are sup- 
posed to exist. It is assumed that the regular 
and irregular components have about the same 
field strength i? « 3 fiG and that both decrease 
exponentially with a scale hight of 1 kpc. The 
scale length of the irregularities is estimated to 
about Lirr ~ 50 pc. 

The magnetic field lines of the loop and fila- 
ment structures are open and perpendicular to 
the galactic disk. In these structures of the disk 
the confinement of cosmic rays is not strong and 
the particles are leaking rapidly. The authors 
assume that the particles leak at this structures 
with diffusive motions along the magnetic field 
lines with a diffusion coefficient Addition- 
ally, outflows with galactic winds with a velocity 
Vg ^ 5 ■ 10~^ pc/a are anticipated. The leakage 
of cosmic rays from the Galaxy is described with 
an one-dimensional advective-diffusion equation 

dn dn d / dn\ ^ 



The diffusion equation is solved numerically 
and the residence time r/j of cosmic-ray parti- 
cles in the galactic disk is obtained. It is found 
that tr depends on the charge of the particles 
and that it is smaller for light nuclei. This im- 
plies that the average mass of cosmic-ray par- 
ticles is expected to increase with energy. The 
relative abundances and the spectral indices for 
each cosmic-ray component are taken from di- 
rect measurements and the all particle spectrum 
is normalized to direct measurements at 10^^ eV. 

The resulting energy spectra for five groups 
of elements are presented in Figure El together 
with the all-particle spectrum obtained. The in- 
dividual components exhibit a relatively smooth 
knee structure, the bends extend over several 
decades in primary energy. The observed icnee 
structure in the all-particle spectrum is not well 
reproduced by this model. The shape of the av- 
erage measured energy spectrum (solid line in 
Figure [T^ is described only moderately by the 
model calculations (filled points). 

Since the heavy component clearly dominates 
above 10^ GeV, a relative strong increase of 
the mean logarithmic mass with energy is ex- 
pected, see Figure jH At the highest energies 
{Eo w 10^ GeV) values of (In^) > 3 are ob- 
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Figure 16: Energy spectra for groups of elements for the model by Roulet and the poly-gonato 
model, see caption of Figure El 



tained. The increase occurs with an almost con- 
stant slope in the logarithmic plot, while the two 
other diffusion models considered, exhibit a dis- 
tinct change in the ioiee region. 

3.12 Diffusion and drift 

Similar to the models discussed previously, 
Roulet et al. consider the drift and diffu- 
sion of cosmic-ray particles in the regular and 
irregular components of the galactic magnetic 
field. Again, a three-component structure of the 
magnetic field is assumed. The regular compo- 
nent is aligned with the spiral arms, reversing 
its directions between consecutive arms. This 
field (with strength Bq) will cause particles with 
charge Z to describe helical trajectories with a 
Larmor radius Rl = p/ (ZeBo). Secondly, a ran- 
dom component is assumed. This will lead to a 
random walk and diffusion along the magnetic 
field direction, characterized by a diffusion co- 
efficient D|| cx i?™. The diffusion orthogonal 
to the regular magnetic field is typically much 
slower, however the energy dependence of Dj_ 
is similar to -Dy. The third component is the 
antisymmetric or Hall diffusion, which is associ- 
ated with the drift of cosmic rays moving across 
the regular magnetic field. The antisymmetric 



diffusion coefficient is Da ~ tlc/3 oc E. 

In this model, for energies below Z-E^ perpen- 
dicular diffusion is the dominant effect, while for 
larger energies the drifts are responsible for the 
dominant escape mechanism. The transition be- 
tween these two regimes is naturally understood 
from the different energy dependence of the two 
diffusion coefficients D± and Da- Both coeffi- 
cients are comparable, Dj_ ~ Da, at an energy 
of the order of Z ■ E/^. This generates a break in 
the all-particle spectrum. 

For the calculations a source spectrum 
dQ/dE oc E~^% with a constant index for all 
species ag = 2.3 is assumed. Below the knee, 
where transverse turbulent diffusion dominates, 
the observed spectral index will be a « 0^ + 1/3, 
while in the drift dominated region above E-Ej., 
a ~ -|- 1 is obtained. This results in a rigidity 
dependent cut-off for the individual elements, as 
can be seen in Figure 1161 For the calculations, 
the cosmic-ray composition has been normal- 
ized to direct measurements at energies below 
10^^ eV. The change of the slope for individual 
components Aa ~ 2/3 is softer than compared 
to the poly-gonato model, where A7 = 2.10 has 
been obtained. 

Nonetheless, the overall agreement between 



21 




i\ I n 1 1 1 1'[ P I I I ui 



5 6 7 8 9 

10 10 10 10 10 

Energy [GeV] 

Figure 17: Energy spectra for groups of elements for the model by Tkaczyk |3Uj (symbols) and the 
poly-gonato model, see caption of Figure El 



the two approaches is quite good. The iron com- 
ponent is clearly dominant at the highest ener- 
gies shown. This results in a relatively strong 
increase of the mean logarithmic mass as func- 
tion of energy, as can be seen in Figure |3 At 
10^ GeV (In^) « 3.75 is reached. 

3.13 Photo-disintegration and diffu- 
sion 

Interactions of cosmic rays with various back- 
ground particles are considered in the next three 
models described. 

Following an idea suggested by Hillas [IH] . 
several authors discuss the possibility of photo- 
disintegration of nuclei in a dense field of soft 
photons around discrete galactic sources, e.g. 
pulsars. Cosmic rays are assumed to be mag- 
netically confined near the source and, there- 
fore, accumulate a large column density on their 
pass across the photon field. These models are 
inspired by the indication in some experiments 
that the mass composition becomes lighter be- 
yond the knee. The knee is explained by the 
onset of photodisintegration as well as due to 
leakage from the Galaxy in the diffusion process 
in the galactic magnetic field. 



Tkaczyk [30] assumes the individual spectra 
to be described by power laws with a spectral 
index of 7p = —2.75 for protons and 'jx = —2.55 
for all other nuclei up to iron using the abun- 
dances observed at 100 TeV by the JACEE ex- 
periment. The photon background is taken to 
have a Planck type distribution with /c^T = 
20 eV. Three processes of energy loss are taken 
into account: pair production, pion photo pro- 
duction on nucleons, and photo-disintegration of 
nuclei. The magnetic field is assumed to be dom- 
inated by its turbulent component namely 93% 
of 31.6 /iT and the radius of the halo is extended 
to 15 kpc. Trajectories of particles have been 
calculated starting at random positions inside 
the galactic disk. 

The energy spectra obtained for elemental 
groups are presented in Figure 1171 The en- 
ergies as given in ,30j have been rescaled by 
6e = —12% in order to match the all-particle 
spectrum (filled circles) with the average mea- 
sured values (solid line). The shape of the all- 
particle spectrum obtained is in good agreement 
with the measured values. The cut-off behavior 
for individual element groups is softer than in 
the poly-gonato model. The steepening in the 
spectrum of primary protons at ~ 10^ GeV is 
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Figure 18: All-particle energy spectrum (left-hand) and mean logarithmic mass versus primary energy 
(right-hand) for the model by Candia et al. (symbols), see also captions of Figures |3] and El 

Results for different values for the temperature T and the column density peL of the photon field are 
shown. 



due to their leakage from the Galaxy. In the 
region between 1 PeV and 30 PeV the iaiee in 
the all-particle spectrum is explained as due to 
photo-disintegration of nuclei and due to leak- 
age. The photo-disintegration process delivers a 
large contribution of secondary protons which 
leads protons to become the dominant mass 
group above 10^ GeV. 

The predicted mean logarithmic mass, shown 
in Figure ^ exhibits a slight increase between 
10^ and 10'' GeV in fair agreement with the 
measurements. Due to the strong contribution 
of secondary protons, {In A) decreases above 
lO''' GeV in contrast to the further increase in- 
dicated by the measured values. 

3.14 Photo-disintegration by optical 
and UV photons 

Photodisintegration of cosmic-ray nuclei in 
dense photon fields around compact sources is 
also discussed in the model by Candia, Epele, 
and Roulet |32j . They investigate in detail 
the potential influence of cosmic rays passing 
through columns of optical and soft UV photons 
in the source region. They vary the column den- 
sity and ksT in the Planck spectrum to study 
the effect on the iaiee structure. Increasing the 
temperature T shifts the knee towards lower en- 
ergies whilst increasing the column density peL 



intensifies the steepening at the iaiee. The au- 
thors conclude that a comparatively small col- 
umn density of about 5 • 10^"^ eV/cm^ and a soft 
photon spectrum with ksT ~ 2 eV describes the 
knee structure best and can, at least partly, be 
responsible for it. 

The total energy spectrum as calculated by 
the authors and the appertaining mean logarith- 
mic mass are plotted in Figure 1181 The spec- 
tra at the source are assumed to be featureless 
extrapolations of the spectra below the knee. 
Spectra for two values of the Planck tempera- 
ture ksT = 1.8 eV and ksT = 10 eV as weh 
as two column densities peL = 5 • 10'^'^ eV/cm^ 
and PeL = 10^^ eV/cm^ are presented. The 
main mechanism of energy loss for nuclei with 
energy Eq < 10^^ eV is the process of photo- 
disintegration. If the typical energy of the pho- 
tons is in the optical range (1-10 eV), the photo- 
disintegration of cosmic-ray nuclei with mass 
A will start to be efficient at energies Eq > 
A ■ 10^^ eV. The appearance of the knee is ulti- 
mately a threshold effect. 

Up to 10^ GeV the shape of the all-particle 
spectrum is compatible with the average mea- 
sured flux. Candia et al. remark that excessive 
fluxes at energies above 10^ GeV are not trouble- 
some since leakage from the Galaxy or different 
efficiencies in the acceleration mechanisms have 
not been considered in the present model, but 
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they can not be disregarded. 

The mean logarithmic mass as derived from 
the present model is shown as function of energy 
in Figure [THl on the right-hand panel. Again, re- 
sults taking into account different temperatures 
and column densities are presented. The photo- 
disintegration of heavy nuclei yields a strong de- 
crease of (In^) above the knee. The extreme 
light composition above 10^ GeV is significantly 
below the observed values. 

3.15 Neutrino interactions in the 
galactic halo 

Dova, Epele and Swain [HJ link the origin of 
the knee to a "GZK-cut-off like effect" of cos- 
mic rays interacting with massive neutrinos in 
the galactic halo. The average number density 
rii, = 337/cm^ of standard-model neutrinos with 
mass nil, < 1 MeV as predicted by Big Bang 
cosmology is assumed to be strongly increased 
due to gravitational clustering in galaxies to a 
value rii, = 1.4 • W^/cnr^. In addition, a mag- 
netic dipole moment fi^, = 5.4 it 0.6 • 10~^ j^B^ 
is adopted for massive neutrinos in order to in- 
crease the cross-section for the inelastic scatter- 
ing of nucleons on the neutrino background. 

Pion production with a threshold energy of 
Ep = 3 PeV is anticipated via the processes 
p + I' I' + A and A — > p + vr. This fixes 
the neutrino mass at nii, = 100 eV. The rela- 
tive heavy mass is attributed to muon or tau 
neutrinos, but exotic objects like hypothetical 
dark matter particles are named as candidates 
as well. These particles are proposed to be in- 
side a spheroidal galactic halo with a radius of 
10 kpc. The cosmic-ray spectrum is modeled 
as a two component mixture (w 60% protons + 
40% iron) with a spectral index 7 = —2.8 for 
both components. The propagation is described 
by a diffusion model, taking into account the 
galactic magnetic field and a residence time of 
about 3 • 10^ a, independent of energy. 

The spectra calculated for the two groups 
"proton" and "iron" are given in a small energy 
range from 1 PeV to 30 PeV in the publication 
and are shown in Figure Hill The energy scale of 
the model has been renormalized hy 5e = —6% 

'^fiB- Bohr magneton 
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Figure 19: Energy spectra for groups of elements 
for the model by Dova et al. j^I] (symbols) and 
the poly-gonato model (lines). The numbers in- 
dicate the range of nuclear charge used for the el- 
emental groups, the all-particle spectrum is rep- 
resented by the filled points and the solid line, 
respectively. 

in order to fit the average all-particle flux of the 
measurements (solid line). In the figure the two 
components are compared with the poly-gonato 
model in two ways: the light component is taken 
as the sum of protons and helium on the one 
hand and as the sum for 1 < Z < 8 on the other 
hand. The remaining elements up to Z = 92 
are used for the heavy component in each case. 
Neither set of spectra is compatible with the cal- 
culations of Dova et al. as can be inferred from 
the figure. The flux of light elements above the 
icnee is overestimated by the model. 

The mean logarithmic mass, as calculated 
from the two groups is almost constant, with an 
average value of about 2 and a slight modulation 
around 10^'^^ GeV, see Figure^ At the highest 
energies shown, the trend towards a light mass 
composition is not supported by the experimen- 
tal values. 

3.16 New physics in the atmosphere 

A completely different reason for the icnee is dis- 
cussed by Kazanas and Nicolaidis j^S]- They 
investigate nucleon interactions hitherto unac- 
counted for in air showers, which transfer en- 
ergy into particles not observable (or not yet 
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Figure 20: All-particle energy spectrum accord- 
ing to calculations by Kazanas & Nicolaidis for 
a general model jSSI and for graviton produc- 
tion j34j (symbols). The average measured all- 
particle flux is shown by the solid line and above 
10® GeV as dashed line. 

observed) by air shower experiments. Candi- 
dates for new physics are supersymmetry, with 
the lightest supersymmetric particle carrying 
away the missing energy, technicolor models 
which produce techni-hadrons not observed by 
present experiments, or models in dimensions 
> 4 lowering the characteristic energy for gravi- 
tation to TeV energies and thus producing gravi- 
ton energy carried away by "invisible" particles. 
Hence, the true energy of the shower inducing 
particle is underestimated in air shower mea- 
surements above the threshold energy. Being 
a threshold effect the new interaction entails a 
relatively sharp knee structure. 

The conventional proton-proton total cross- 
section rises slowly with energy. In the model it 
is considered to have a constant value of ctq = 
80 mbarn. The additional cross-section an{E) 
describes the proposed new type of interaction, 
rising strongly with energy. At energies above 
the threshold a fraction of protons interacts with 
the probability P„(E) = an{E)/[(Ja{E)+an{E)] 
through the new channel. Hence, more and more 
energy is missing in the measured components. 

As primary all-particle energy spectrum the 
authors assume a single power law with a 
cut-off at an energy Eq of the form / = 
S-Texp(-S/£;o) with 7 = 2.75 and Eq = 



-|^q18.5 jj-^ addition, at energies above 10^® eV, 
an extragalactic component is introduced. The 
all-particle spectrum above the atmosphere is 
modulated through the proposed new interac- 
tions. The resulting apparent spectrum, re- 
constructed from observations below the atmo- 
sphere is plotted in Figure 1201 (circles). The 
model produces a spectrum with a shape sim- 
ilar to the average measured values, the knee 
occurs at a relative high energy above 10'' GeV. 
The authors point out that the present model 
is very simple. Hence, some deviations between 
measurements and the predicted spectrum are 
possible. 

In a second article Kazanas and Nicolaidis 
perform more detailed calculations of the cross- 
section in the case of graviton radiation ■ 34 . As 
new physics in the interactions low scale gravi- 
tation is discussed. The graviton is considered 
to propagate not only in the usual four dimen- 
sions but in an additional compactifyed large 
dimension 6. Thus the fundamental (Planck) 
scale of gravity Mf is reduced to TeV ener- 
gies. Any scattering process at center-of-mass 
energies E > Mj should be accompanied by 
abundant graviton production, i.e. collinear 
bremsstrahlung of gravitons. The cross-section 
for graviton production in pp collisions is calcu- 
lated to be an oc {^/s / M j)^'^^ . A fit to the mea- 
sured cosmic-ray spectrum yields Mj 8 TeV 
and (5 = 4. This indicates a strong increase of 
the graviton production cross section as func- 
tion of energy. Even at energies just above Mf, 
graviton radiation will be an important process. 

Again, as primary spectrum a single power 
law has been assumed and the apparent spec- 
trum has been calculated. It is shown in Fig- 
ure 1201 (squares). The spectrum is very similar 
to the one obtained from the more general con- 
sideration for different processes as mentioned 
above. 

4 Discussion 

All models discussed in this article are summa- 
rized in Table [21 Most of them yield similar 
all-particle spectra. On the other hand, the in- 
dividual element spectra are quite different in 
some models. Some predict a very weak knee 
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Table 2: Synopsis of all models discussed. 

Author(s) 



Model 

Source /Acceleration: 

Acceleration in SNR 

Acceleration in SNR + radio galaxies 

Acceleration by oblique shocks 

Acceleration in variety of SNR 

Single source model 

Reacceleration in the galactic wind 

Cannonball model 
Propagation/Leakage from Galaxy: 

Minimum pathlength model 

Anomalous diffusion model 

Hall diffusion model 

Diffusion in turbulent magnetic fields 

Diffusion and drift 
Interactions with background particles: 

Diffusion model + photo-disintegration 

Interaction with neutrinos in galactic halo 

Photo-disintegration (optical and UV photons) 
New interactions in the atmosphere: 

Gravitons, SUSY, technicolor 



structure, like e.g. the acceleration in cannon 
balls by Plaga and the anomalous diffusion de- 
scribed by Lagutin et al, while others obtain a 
steep falling flux at the icnee, an example is the 
diffusion model by Kalmykov and Pavlov. The 
results from indirect measurements seem to fa- 
vor a relative steep cut-off, see section |2l 

Most models predict featureless spectra with 
a simple shape, however, some scenarios yield 
distinct features, which may have different the- 
oretical origins. Kalmykov and Pavlov as well as 
Swordy consider a diffusion model and a Leaky 
Box model, respectively and obtain bumps at 
the individual knee positions. Consequently, 
they find a dip at iaiee energies in the mean 
logarithmic mass. The approaches by Svesh- 
nikova (acceleration in supernova remnants) as 
well as Volk and Zirakashvili (reacceleration in 
the galactic halo) lead both to a twofold shape 
of the element spectra. Yet, also common fea- 
tures can be recognized, like the relative sharp 
kinks in the spectra as found by Erlykin and 
Wolfendale as well as by Stanev et al., both mod- 
els describe acceleration in supernova remnants. 
The distinct shape leads in both cases to struc- 



Berezhko &: Ksenofontov ^S] 
Stanev et al. [T9] 
Kobayakawa et al. |2Uj 
Sveshnikova |2j 
Erlykin & Wolfendale [H] 
Volk and Zirakashvili [22] 
Plaga mi 

Swordy 

Lagutin et al. [2^1 

Ptuskin et al. [211, Kalmykov & Pavlov [121 
Ogio and Kakimoto [2H1 
Roulet et al ^ 

Tkaczyk [HOI 
Dova et al. [3J 
Candia et al. [S21 

Kazanas & Nicolaidis |33| l3lj 



tures in the (In^) distribution. 

The review has shown further that even simi- 
lar conceptions lead to (In^) values covering the 
whole range depicted in Figure 0J As an exam- 
ple, acceleration in supernova remnants yields 
{In A) values at the lower edge of the experi- 
mental values (Berezhko and Ksenofontov) and 
also at the upper border (Kobayakawa et al.). 

Unfortunately, the experiments give no con- 
clusive energy spectra for elemental groups (Fig- 
ures n - [31) and the mean logarithmic mass 
covers a relatively wide range, see Figure [Ij 
Most astrophysical informations available from 
air shower measurements depend on the model 
used to describe the interactions in the atmo- 
sphere. Hence, an interpretation and evaluation 
of the data should be treated with care and it 
is very challenging to characterize in a quanti- 
tative way the best model to describe the origin 
of the icnee. 

In models which introduce new type of in- 
teractions in the atmosphere, like the one just 
discussed by Kazanas and Nicolaidis or the ap- 
proaches by Nikolsky et al. |;44^ or by Petrukhin 
|45j . the threshold for a new type of interaction 
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depends on the energy per nucleon. Thus, such 
models ultimately arrive at a mass dependent 
cut-ofF for individual element spectra. The in- 
vestigations of the K AS CADE group however 
seem to indicate a rigidity dependent cut-off for 
individual elements [H 01 IH El El • Also with the 
phenomenological approach of the poly-gonato 
model a rigidity dependent approach is favored 

HI- 

In addition, the simultaneous observations 
of a icnee in the electromagnetic, muonic and 
hadronic component by the KASCADE group 
^, 46 , 47j yield no major inconsistencies between 
the different air shower components. Therefore, 
a proposed new interaction would have to " hide" 
the energy in such a way, that no inconsisten- 
cies occur, when the measurements of different 
shower components are interpreted with stan- 
dard particle physics. From the present exper- 
imental results new particle physics for high- 
energy interactions in the atmosphere seems not 
to be a likely explanation for the icnee. 

The interaction of cosmic rays with back- 
ground particles like neutrinos or photons dur- 
ing their propagation causes fragmentation of 
heavy nuclei and the production of a large 
amount of secondary light nuclei, preferable pro- 
tons. Models based on such an approach, like 
the models by Tkaczyk, Dova et al.^, and Can- 
dia et al. yield a very light mean logarithmic 
mass above icnee energies. On the other hand, 
the measurements indicate an increase of (In A) 
with energy, see Figure^ thus disfavoring inter- 
actions with background particles in the Galaxy 
as an explanation for the icnee. In addition, the 
cross-sections for such interactions depend on 
the energy per nucleon. This yields a mass de- 
pendent cut-off for individual element spectra, 
which seems not to be favored by the measure- 
ments as mentioned above. 

Numerous theories consider diffusive propa- 
gation of cosmic rays through the interstellar 
medium and leakage from the Galaxy as an ex- 
planation of the icnee. The escape probability 
from the Galaxy in these models depends on the 

*Dova et al give their results only for a very limited 
range in primary energy, but above lO'^ GeV a trend to- 
wards an extremely light mass composition can be rec- 
ognized. 



nuclear charge Z and a rigidity dependent cut- 
off behavior is obtained for the individual ele- 
ment spectra. These models yield a more or less 
strong increase of the mean logarithmic mass 
above the icnee as consequence of the leakage 
starting with light nuclei. The change of the 
spectral slope at the icnee is not very distinct 
in the Leaky Box model by Swordy as well as 
in the anomalous diffusion model by Lagutin et 
al. Consequently, only a modest increase of the 
mean logarithmic mass is obtained for energies 
above the icnee, at the highest energies the pre- 
dicted (In A) values are at the lower edge of the 
measurements. 

Propagation of cosmic rays in the Galaxy, tak- 
ing into account contributions from the regular 
magnetic field, a random component and anti- 
symmetric diffusion are described by Kalmykov 
and Pavlov, Ogio and Kakimoto, as well as 
Roulet et al. Although they consider simi- 
lar scenarios, quite different energy spectra are 
obtained. Ogio and Kakimoto derive a very 
smooth transition in the icnee region, while 
Kalmykov and Pavlov as well as Roulet et al. 
calculate a cut-off behavior similar to the mea- 
surements. The mean logarithmic mass ob- 
tained for the three models covers a wide range, 
and around 10^ GeV a significantly different be- 
havior for the individual models is found. On 
the other hand, they all result in similar values 
of (In A) ~ 3 around 10^'^ GeV, well compatible 
with the measurements. 

Several models attribute the icnee to the ac- 
celeration process of cosmic rays. In the can- 
nonball model described by Plaga, the spectral 
slopes change only by about A7 ~ 0.3 at the 
icnee and protons are the dominant nuclei even 
above the icnee. Hence, the resulting mean log- 
arithmic mass is almost a constant function of 
energy. It is the weakest energy dependence of 
all models discussed. In the single source model 
by Erlykin and Wolfendale heavy nuclei domi- 
nate above the icnee and the {In A) values are 
in fair agreement with the measurements. This 
can also be stated for the concept of Volk and Zi- 
rakashvili 1221, which allows for (re) acceleration 
of cosmic rays up to energies of the ankle. 

Diffusive shock acceleration in supernova rem- 
nants is considered by Berezhko and Kseno- 
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fontov, Stanev et al., Kobayakawa et al., and 
Sveshnikova. Different variations of the basic 
theory of Fermi acceleration have been appUed. 
Due to the finite hfetime of the shock front 
the maximum energy reached is Umited. Typ- 
ical values of Z ■ (0.1 — 3) PeV are attained. 
Kobayakawa et al. could increase the maxi- 
mum energy reached in their calculations by 
the investigation of special magnetic field con- 
figurations. Sveshnikova takes into account re- 
cent supernova observations, which indicate the 
existence of very energetic explosions (hyper- 
novae). In their high- velocity shock fronts par- 
ticles can effectively be accelerated to high en- 
ergies. The spectra obtained by the latter two 
scenarios are similar to the measurements and 
the poly-gonato model. 

The poly-gonato model is purely phenomeno- 
logical. The spectra obtained are most com- 
patible with the acceleration in supernova rem- 
nants as described by Sveshnikova, especially 
the heavy enriched case, and diffusive prop- 
agation through the Galaxy as discussed by 
Kalmykov and Pavlov as well as Roulet et al. 
Hence, these models could be a good theoretical 
foundation of the phenomenological ansatz. 

Most models discussed represent principle 
ideas, dealing with individual aspects of cosmic- 
ray acceleration, propagation, or interactions. 
Most probably, combinations of this basic pro- 
cesses are needed to describe nature realistically 
and more precise data will be necessary for a 
fine tuning of the free parameters in the models. 
These parameters have been adjusted to differ- 
ent sets of experimental results in the individual 
models. Most likely some of their predictions 
change, if a different set of experimental results 
is used. 

At present, the experimental results do not al- 
low to select a single approach as the best model 
to describe nature. The origin of the knee in 
the all-particle energy spectrum is not yet ex- 
plained. However, within the boundaries given 
by the experiments, acceleration in supernova 
remnants and diffusive propagation through the 
Galaxy seem to be very attractive models to un- 
derstand the origin of the knee. 
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